KU-PH-013 
TTP 12-042 
SFB/CPP-12-85 

Full 0(a) electroweak radiative corrections to 
e + e~ -> ttj with GRACE-Loop 



P.H. Khiem" 4 '^, J. Fujimoto" 4 , T. Ishikawa A , T. Kaneko" 4 , 

K. Kato c , Y. Kurihara A , 
Y. Shimizu" 4 , T. Ueda D , J. A.M. Vermaseren^, Y. Yasui F 

A ^KEK, Oho 1-1, Tsukuba, Ibaraki 305-0801, Japan. 
B > SOKENDAI University, Shonan Village, Hayama, Kanagawa 240-0193 Japan. 
c ^ Kogakuin University, Shinjuku, Tokyo 163-8677, Japan. 
D ^ Karlsruhe Institute of Technology (KIT), D-76128 Karlsruhe, Germany. 
E ^NIKHEF, Science Park 105, 1098 XG Amsterdam, The Netherlands. 
^ Tokyo Management College, Ichikawa, Chiba 272-0001, Japan. 



Abstract 

We present the full 0(a) electroweak radiative corrections to the process e + e~ —> 

tt^y at the International Linear Collider (ILC). The computation is performed with 

the help of the GRACE-Loop system. We present the total cross-section and the 

top quark forward-backward asymmetry (Aps) as a function of the center-of-mass 

energy and compare them with the process e + e~ —> tt. We find that the value 

of Afb in tt'y production is larger than Aps in ti production. It is an important 

result for the measurement of the top quark forward-backward asymmetry at the 

ILC. Applying a structure function method, we also subtract the QED correction 

to gain the genuine weak correction in both the a scheme and the scheme (6yy). 

G 

We obtain numerical values for 5-^ which are changing from 2% to —24% when we 
vary the center-of-mass energy from 360 GeV to 1 TeV. 



1 Introduction 



The experimental results of CDF pQ and DO [2] on the measurement of top pair produc- 
tion at the Tevatron show an unexpected large top quark forward-backward asymmetry. 
The precise theoretical calculations of the top pair production play an important role 
in explaining the experimental data. QCD radiative corrections to top pair production 
from proton-proton collisions were calculated by several authors [3J, [Jj. However, the 
measurement is affected by a huge background from QCD. A good example is the gg — > ti 
reaction. In the future, the measurement will be performed at the ILC without QCD 
background. Therefore, we consider the precise calculations of top pair production and 
top pair with photon production in e + e~ collisions. A completed full one-loop electroweak 
correction calculation to the process e + e~ — > ti has already been presented in refs [3], [5J, 
[7j. In this paper, we calculate the full 0(a) electroweak radiative corrections to both the 
process e + e~ — > ti and e + e~ — > tt'y at the ILC. 

The data of the ATLAS [8] and CMS [£] experiments prove the existence of a new 
boson with mass around 126GeV. It is assumed to be the standard-model Higgs particle. 
Once the discovery of the Higgs boson is confirmed, the next important task is to measure 
its properties. However, it is clear that such a measurement is much easier at the cleaner 
environment of the ILC than at the LHC with its large QCD backgrounds. To measure the 
properties of the new boson it is important that also the radiative correction calculations 
for the ILC take the complete standard model into account. 

The experiments at the ILC require a precise determination of the luminosity which 
will be based on higher order theoretical calculations of the Bhabha scattering cross- 
section. Thus, the computation of electroweak radiative corrections to the process e + e~ — > 
e + e _ 7 is mandatory. This will be our eventual target. However, as a first step, we are 
going to calculate the process e + e~ — > tt^y which is easier in several respects: there are 
fewer diagrams and the numerical cancellations between the diagrams are less severe. It 
will provide a framework for our target calculation. 

In the scope of this paper, we discuss the full 0(a) electroweak radiative corrections 
to the process e + e~ — > ti'-f at ILC. We then examine the numerical results of the top 
quark forward-backward asymmetry as well as the genuine weak corrections in both the 
a scheme and the scheme [TO] as compared to the process e + e~ — > ti. 

The paper is organized as follows. In section 2 we introduce the GRACE-Loop system 
and set up the calculation. In section 3 we discuss the numerical results of the calculation. 
Future plans and conclusions of our paper are presented in section 4. 
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2 GRACE-Loop and the process e e + — > tt^y 
2.1 GRACE-Loop 

The computation is performed with the help of the GRACE-Loop system which is a generic 
program for the automatic calculation of scattering processes in High Energy Physics. Of 
course, with a system this complicated still under development, it is important to have as 
many tests as possible of the correctness of the answer. Hence the GRACE-Loop system 
has been equipped with non-linear gauge fixing terms in the Lagrangian which will be 
described in some of the next paragraphs. The renormalization has been carried out with 
the on-shell renormalization condition of the Kyoto scheme, as described in ref [11]. The 
program was presented and checked carefully with a variety of 2 — > 2-body electroweak 
processes in ref [12J. The GRACE-Loop system has also been used to calculate 2 — > 3- 
body processes such as e + e~ — )■ ZHH [13], e + e~ tiH [H], e + e~ — ?• vvH [15] and the 
2 — > 4-body process e + e~ — > jiuud [16], e + e~ — > v^v^HH [17] . 

The steps of calculating a process in the GRACE system are as follows. First the 
system requires input files that describe the Feynman rules of the model. In this case, we 
use the standard model. These files are considered part of the system but for different 
models the user would have to provide them. Next a (small) file is needed that selects 
the model, the names of the incoming and outgoing particles, and one of a set of prede- 
fined kinematic configurations for the phase space integration. In the intermediate stage 
symbolic manipulation handles all Dirac and tensor algebra in n- dimensions, reduces the 
formulas to coefficients of tensor one-loop integrals and writes the formulas in terms of 
FORTRAN subroutines on a diagram by diagram basis. For this manipulation either 
FORM [H] or REDUCE [H] is used. The FORTRAN routines will be combined with 
libraries which contain the routines that reduce the tensor one-loop integrals into scalar 
one-loop integrals. The scalar one-loop integrals will be numerically evaluated by one 
of the FF [2U] or Loop Tools [2TJ packages. The ultraviolet divergences (UV-divergences) 
are regulated by dimensional regularization and the infrared divergences (IR-divergences) 
will be regulated by giving the photon an infinitesimal mass A. Eventually all FORTRAN 
routines are compiled and linked with the GRACE libraries which include the kinematic 
libraries and the Monte Carlo integration program BASES [22] . The resulting executable 
program can then calculate cross-sections and generate events. 

Ref [12] describes the method used by the GRACE-Loop system to reduce the tensor 
one-loop five- and six-point functions into one-loop four-point functions. 

The GRACE-Loop system allows the use of non-linear gauge fixing conditions [23] 
which are defined by: 

C GF = -— - leaA^ - igc w pZ^)W» + + Z w 9 -(v + ~5H + iR X z)x + \ 2 
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- — {d ■ Z + £ * (« + sH) X3 f - —{d ■ Af . (1) 

We are working in the 't Hooft-Feynman gauge by choosing ^ w — £ z = £4 = 1. Therefore, 
there is no contribution of the longitudinal term in the gauge propagator. This choice has 
not only the advantage of making the expressions much simpler. It also avoids unneces- 
sary large cancellations, high tensor ranks in the one-loop integrals and extra powers of 
momenta in the denominators which cannot be handled by the FF package. 

The GRACE-Loop system can also use an axial gauge for external photons. This has 
two advantages. 

1. It cures a problem with large numerical cancellations. This is very useful when 
calculating the process at small angle and energy cuts of the final state particles. 

2. It provides a useful tool to check the consistency of the results which, due to the 
Ward identities, must be independent of the choice of the gauge. 

2.2 The numerical test of the process e e + — > ti^y 

The full set of Feynman diagrams with the non-linear gauge fixing as described before 
consists of 16 tree diagrams and 1704 one- loop diagrams (of which 168 are pentagon 
diagrams). In Fig [T] we show some selected diagrams. 

The results are checked carefully by three kinds of consistency tests. These tests are 
performed with quadruple precision at a few random points in the phase space. The 
first test is ultraviolet finiteness of the results. This test is done on all virtual one-loop 
diagrams and their counter terms and we treat Cuv — l/ e ~ 7_B + log47r as a parameter. In 
order to regularize the infrared divergences, we give the virtual photon a fictitious mass, 
A = 10~ 17 GeV. In the tabled] we present the numerical results of the test at one random 
point in the phase space. The result is stable over more than 30 digits for various values 
of the ultraviolet parameter. 



Cuv 




23ft(7jy. ee 7L 00 p) 







-5.3131630^ 


35402176811947711662^ 


33 1 7605E- 3 


10 


-5.3131630^ 


35402176811947711662^ 


33 1 7726E" 3 


100 


-5.3131630^ 


3540217681 1947711662^ 


332 4 04£- 3 



Table V. Test of Cuv independence of the amplitude. In this table, we take the non-linear 
gauge parameters to be 0, A = 10 _17 GeV and we use 1 TeV for the center-of-mass energy. 

The second test is the independence of the result on the fictitious photon mass A. In 
this case, we take Cuv — 0- This test will be performed by including as well the virtual 
loop diagrams as the soft bremsstrahlung contribution. In table [2] the numerical results 
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Graph 933 Graph 1159 Graph 1291 




produced by GRACEFICJ 



Figure 1: Typical Feynman diagrams as generated by the GRACE-Loop system. 

of the test are presented. We find that the result is stable over more than 15 digits when 
varying the parameter A over a wide range. 

The independence of the result on the five parameters a, (3, 5, R, e is also checked. The 
result is presented in table [3J We find that the result is stable over more than 26 digits 
while varying the non-linear gauge parameters. 

Finally, we check the stability of the result versus the soft photon cut parameter (k c ). 
This test includes both the soft photon and the hard photon contributions. The hard 
photon bremsstrahlung part is the process e + e~ — > ttjj. It is important to note that 
we have two photons at the final state. One of them has to be applied an energy cut of 
E° ut > 10 GeV and an angle cut of 10° < 9° ut < 170°. Another one is a hard photon 
with energy is greater than k c and smaller than the first photon's energy. This part will 
be generated by the tree level version of GRACE [21] with the phase space integration by 
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A [GeV] 


23?(72Y ee 7L 00p )+soft contribution 


io- 17 


-1.6743892369492021873805611201763810£- 3 


lO" 19 


-1.6743892369492020397654354220438766E- 3 


IO" 21 


-1.6743892369492020382892402083349623E- 3 



Table 2: Test of the IR finiteness of the amplitude. In this table we take the non-linear 
gauge parameters to be 0, Cuv — and the center-of-mass energy is 1 TeV. 



(a, f3, k, 6, e) 


23? ( 7jy ee TLoop ) 




(0,0,0,0,0) 


-5.3131630854021768119477116628317605-ET 


-3 


(1,2,3,4,5) 


-5.3131630854021768119477116637537265^ 


-3 


(10,20,30,40,50) 


-5.3131630854021768119477116582762373^ 


-3 



Table 3: Gauge invariance of the amplitude. In this table, we set Cuv = 0, the photon 
mass is 10 _17 GeV and a 1 TeV center-of-mass energy. 

BASES. The result is tested by changing the value of k c from 10- 5 GeV to 0.1 GeV. In 
table H] we find that the results are in agreement with an accuracy which is better than 
0.1% when we vary k c . 



k c [GeV] 






0"S 




&S+H 


IO" 5 


4.172723£- 


-02 


5.885469£~ 


-02 


0.10058192 


IO" 3 


2.926684£~ 


-02 


7.131737.B- 


-02 


0.10058421 


10- 1 


1.678994E" 


-02 


8.377319£- 


-02 


0.10056313 



Table 4: Test of the fc c -stability of the result. We choose the photon mass to be 10 
GeV and the center-of-mass energy is 1 TeV. The second column presents the hard photon 
cross-section and the third column presents the soft photon cross-section. The final column 
is the sum of both. 

We found that the numerical results are in good agreement when varying Cuv, the 
gauge parameters, photon mass, and k c . Hereafter, we set A = IO -17 GeV, Cuv — and 
a = (3 = 5 = k = e = 0. 

3 Results 

Our input parameters for the calculation are as follows. The fine structure constant in the 
Thomson limit is a -1 = 137.0359895. The mass of the Z boson is Mz = 91.187 GeV. In 
the on-shell renormalization scheme we take the mass of the W boson (M^y) as an input 
parameter. It will be derived through the electroweak radiative corrections to the muon 
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decay width (Ar) [2SJ with = 1.16639 x 1(T 5 GeV~ 2 . Therefore, M w is a function 
of Mjj- In this calculation, we take Mh = 120 GeV and the numerical value of M\y is 
80.3759 GeV. For the lepton masses we take m e = 0.51099891 MeV, m T = 1776.82 MeV 
and m M = 105.658367 MeV. For the quark masses we take m u = 1.7 MeV, rrid = 4.1 
MeV, m c = 1.27 GeV, m s = 101 MeV, m t = 172.0 GeV and m b = 4.19 GeV. We apply 



an energy cut of E° u > 10 GeV and an angle cut of 10° < 9™ < 170° on the photon. 

All numerical results are generated by the GRACE-Loop system. For ti production 
the results were first checked with the results in refs [5] , [6] , [7] . Then we use the values of 
the parameters above to produce the results of ti production in this paper and compare 
them with tiy production. 

In Fig [2] the total cross-section is a function of the center-of-mass energy ^/i. We vary 
the value of y/s from 360 GeV to 1 TeV. We find that the cross-section is largest near 
the threshold, y/s around 550 GeV for tij production and 410 GeV for ti production. 
The total cross-section of tiy production is considerably less than 10% of the total cross- 
section for the ti reaction. In addition we find a negative correction for tij production in 
contrast to the positive correction for ti production. 
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Figure 2: The total cross-section as a function of center-of-mass energy. The left figure is 
the result of ti production and the right figure shows the result of the tij reaction. The 
triangle points are the result of the tree level calculation while the rectangular points are 
the sum of the tree level calculation combined with the full one-loop electroweak radiative 
corrections. Lines are only guide for the eyes. 
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The full 0{a) electroweak corrections take into account the tree graphs and the full 
one-loop virtual corrections as well as the soft and hard bremsstrahlung contributions. 
The relative correction is defined as 

oew = 1- [*) 

&Tree 

In order to extract the genuine weak correction in the scheme, we first evaluate 
the QED initial radiative correction (Sqed)- Applying the structure function method 
described in ref Sqed is defined as 



r &QED ~ °T ree 

OQED = ; 

CTree 

with 



(3) 



1 

Vqed(s) = J dx K(x, s) cr (s(l - x)), (4) 



here H(x, s) is a radiator which is defined by formula (11.213) in ref [26] : 

H(x,s) = Afix?- 1 - (3(1 - |) 



(3 2 r l + 3(l-a^ 2 



+ — 



4(2-x)lnx — ln(l - x) - 6 + x (5) 



x 



with (3 = & (ln(^) - l) and A = 1 + a (fln(^) + £ - 2). 

After obtaining the QED correction, we define the genuine weak correction in the a 
scheme: 

5w = Sew — Sqed- (6) 

Having subtracted the genuine weak correction in the a scheme, one can express the 
correction in the scheme. Next we subtract the universal weak correction which is 
obtained from Ar. The genuine weak correction in the G M scheme is defined by 



tify =5 W - nAr, (7) 

with Ar = 2.55% for M# = 120 GeV and n = 3(2) for ttj (for tt) production respectively. 

In Fig [21 we present the full electroweak correction and the genuine weak correction in 
both the a and the G^ schemes for ti'y production as compared to tt production. These 
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corrections are shown as a function of the center-of-mass energy, y/s. We vary y/s from 
360 GeV to 1 TeV. The figures show clearly that the QED correction is dominant in the 
low energy region. In the high energy region it is much smaller (~ —5% at 1 TeV). In 
contrast to the QED correction the weak correction in the a scheme is less than 10% for 
low energies but reaches —16% at 1 TeV center-of-mass energy. For tt^f production, we 
find that the value of the genuine weak correction in the scheme varies from 2% to 
-24% over yfi from 360 GeV to 1 TeV. 



e e + -)-tt e e + ->■ tt^y 




300 400 500 600 700 800 900 1000 300 400 500 600 700 800 900 1000 



Center-of-mass energy [GeV] Center-of-mass energy [GeV] 



Figure 3: The full electroweak correction and the genuine weak correction as a function 
of the center-of-mass energy. The left figure shows the results for tt production while the 
right figure shows the results for tt'j production. The circle points represent the QED 
correction, the empty rectangle points are the results for the full electroweak correction 
while the triangle points are the results for the genuine weak correction in the a scheme. 
The filled rectangle points represent the results of the genuine weak correction in the 
scheme. Lines are only guide for the eyes. 

Now we turn our attention to the forward-backward asymmetry A FB . This quantity 
is defined as 

a(0° <6 t < 90°) - (j(90 o < 9 t < 180°) 
FB ~ a(0° <6 t < 90°) + a(90° <6 t < 180°) ' ^ ' 

with 6 t the angle of the top quark. 
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Fig H] shows the results for Aps as a function of the center-of-mass energy. The 
figures show clearly that the top quark asymmetry in the full results is smaller than the 
asymmetry at the tree level results only. 




Figure 4: The top quark forward-backward asymmetry as a function of the center-of- 
mass energy. Left figure is the results for ti production and right one is the results for ti~f 
production. The triangle points represent the tree level results and the rectangle points 
are the results including the full radiative corrections. Lines are only guide for the eyes. 

In Fig [5] we compare the values of Aps i n ti^y production directly with its value for ti 
production. From the figures, we find that Afb m itl production is larger than Apb in ti 
production. This is the most important result of the paper. The effect should be clearly 
observable at the ILC. 



4 Conclusions 

We have presented the full 0(a) electroweak radiative corrections to the process e + e~ — > 
tt'y and e + e~ — > ti at ILC. The calculations were done with the help of the GRACE-Loop 
system. 

The results have been checked carefully by three kinds of consistency tests. The 
results are numerically stable when quadruple precision is used. This shows the power of 
the system. 
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Figure 5: The value of the top quark asymmetry in tt'y production as compared to tt 
production. The rectangular points (circle points) represent the result for tt'y production 
(tt production) respectively. Lines are only guide for the eyes. 

We find that the numerical value of the genuine weak corrections in scheme varies 
from 2% to —24% in the range of center-of-mass energy from 360 GeV to ITeV. We also 
obtain a large value for the top quark forward-backward asymmetry in the tt'y process as 
compared with the one in tt production. 

We also introduce the axial gauge for the external photon in the GRACE-Loop sys- 
tem. It helps to avoid a large numerical cancellation problem. This is very useful when 
calculating Bhabha scattering at small angle and energy cuts of the final state particles. 
Bhabha scattering and related processes are not only used as luminosity monitor, but 
also play an important role as backgrounds for the process e _ e + — > X~X + 7, which is a 
very interesting reaction for the search for dark matter. We will address it in a future 
publication. 

In addition, this calculation will provide a framework for calculating the full O(o) 
electroweak radiative corrections for the process e + e _ — >■ e + e~7. This reaction and these 
corrections will play an important role at future e + e~ colliders like the ILC. 
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